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Abstract 

Cygnus X is a conspicuous massive-star forming region in the Local Arm of the Galaxy at ~ I A kpc from the solar 
system, y-ray observations can be used to trace cosmic rays (CRs) interacting with the ambient interstellar gas and 
low-energy radiation fields. Using the Fermi Large Area Telescope (LAT) we have discovered the presence of a 50-pc 
wide cocoon of freshly-accelerated CRs in the region bounded by the ionization fronts from the young stellar clusters. 
On the other hand, the LAT data show that the CR population averaged over the whole Cygnus complex on a scale 
of ~ 400 pc is similar to that found in the interstellar space near the Sun. These results confirm the long-standing 
hypothesis that massive star-forming regions host CR factories and shed a new light on the early phases of CR life in 
such a turbulent environment. 

Keywords: cosmic rays, gamma rays, acceleration of particles, open clusters and associations 
PACS: 98.20.Af, 98.70.Rz, 98.70.Sa 



1. Introduction 

Cosmic rays (CRs) are one of the most intriguing 
puzzles of the last century Physics. It is strongly ad- 
vocated that Galactic CRs are accelerated by supernova 
remnants (SNRs), which are energetic and numerous 
enough to sustain the CR population directly measured 
near the Earth |1J. A plausible mechanism for CR ac- 
celeration in SNRs is provided by non-linear diffusive 
shock acceleration. The presence of high-energy elec- 
trons in SNRs is undoubtedly demonstrated by their 
multiwavelength spectrum. On the other hand, the ac- 
celeration of nuclei is an elusive phenomenon, which, 
from the observational point of view, so far is supported 
mainly by the impact of efficient ion acceleration on the 
thermohydrodynamics of the Shockwaves, as inferred 
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from X-ray observations, and from the spectral shape 
of y-ray emission [e.g.|2l. 

The isotopic abundances measured for CRs indicate 
that ~ 20% of the material is synthesized by Wolf- 
Rayet (WR) stars 100 kyr before the acceleration jS). 
WR stars represent an evolutionary stage of massive OB 
stars, and the latter cluster in space and time within 
the so-called OB associations. Additionally, ~ 80% 
of SNRs originate from the gravitational collapse of a 
massive-star core, therefore often inside an OB associ- 
ation. All these considerations have given credence to 
the hypothesis that at least part of the CRs are acceler- 
ated by the repeated action of Shockwaves from massive 
stellar winds and SNRs inside massive stellar clusters. 

The propagation of high-energy CRs in the interstel- 
lar space is often described in terms of diffusion by mag- 
netohydrodynamic turbulence with the possible inclu- 
sion of convection and reacceleration Q. While prop- 
agating, CRs loose their energy through many different 
interaction processes. Electromagnetic radiation aris- 
ing from such interactions is the only probe we have 
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of CR properties in the Galaxy beyond direct measure- 
ments performed in the solar system. Interstellar y-ray 
emission is produced by CR interactions with the in- 
terstellar medium, via nucleon-nucleon inelastic colli- 
sions and electron Bremsstrahlung, and with the low- 
energy interstellar radiation fields via inverse-Compton 
(IC) scattering by electrons. 

The intermediate steps, i.e. the escape of CRs from 
their sources and the early propagation in the surround- 
ing medium, have escaped observations so far. Massive 
stars generate intense radiation, powerful winds and ex- 
plosions which alter their neighborhood and power fast 
Shockwaves and strong supersonic turbulence. If part 
of the CRs are accelerated in regions of massive-star 
formation, this turbulent medium influences their evolu- 
tion, e.g. through confinement and reacceleration phe- 
nomena counter-balanced by increased radiative losses. 

The Cygnus X region is home to numerous mas- 
sive stellar clusters Q embedded in a complex of giant 
molecular clouds. Most of the molecular gas along the 
line of sight shows interactions with the stellar clusters, 
notably with Cyg OB2 at a distance of ~ 1.4 kpc ||6|. 
At a comparable distance we find the SNR y Cygni, e.g. 
Q, which appears as a potential CR accelerator In ad- 
dition, Milagro detected very high-energy y-ray emis- 
sion from the complex f8l, pointing to the presence of 
enhanced CR densities. 

In this paper we summarize the results of our anal- 
ysis of two years of high-energy y-ray observations of 
the Cygnus complex by the Fermi Large Area Telescope 
(LAT) S. The results, discussed in detail in IfTOlfTTH . 
reveal an unexpected scenario, where the low-density 
cavities blown by the massive stellar clusters form a co- 
coon for freshly-accelerated particles, whereas the CR 
population averaged over the whole complex is similar 
to that found in other more quiet segments of the Local 
Arm. 

2. Data and Analysis 

We analyzed the region at Galactic longitudes 72° < 
/ < 88° and latitudes -\5°b+\5°. The multiwavelength 
data and analysis method are detailed in |11|: in this 
section we summarize the main features. 

We have used y-ray data from the first two years of 
observations by the Fermi LAT between 100 MeV and 
100 GeV, applying tight rejection criteria to reduce the 
backgrounds due to charged CRs and the Earth's atmo- 
spheric emission. Below 1 GeV we accepted only y- 
rays that converted into pairs in the thin section of the 
LAT tracker, to reduce the angular resolution degrada- 
tion due to multiple Coulomb scattering. 



The model described hereafter was fitted to the LAT 
data by a binned maximum likelihood procedure on a 
0.125° X 0.125° angular grid and independently for 10 
energy bins spanning the whole energy range consid- 
ered, using post-launch LAT instrument response func- 
tions. 

Since the densities of the bulk of Galactic CRs are 
expected to vary only on scales larger than those of the 
gas clouds, we modeled interstellar y-ray emission from 
their interactions as a linear combination of components 
traced by atomic hydrogen (Hi) and CO (as a surro- 
gate of molecular hydrogen, H2) split along the line the 
sight to separate the Cygnus complex from the more dis- 
tant Perseus and outer spiral arms, plus visual extinction 
residuals to take into account the neutral interstellar gas 
missed by these two tracers. We added an IC model cal- 
culated using the GALPROP code f4| and an isotropic 
component to account for the extragalactic diffuse y-ray 
background and the residual background due to charged 
CR interactions in the LAT. 

Below a few GeV the y-ray emission from the region 
is dominated by the three bright pulsars J2021 -1- 3651, 
J2021 + 4026 and J2032 H- 4127. To be more sensi- 
tive to interstellar structures, at energies < 10 GeV their 
emission was dimmed by excluding from the analysis 
photons detected in their periodic pulses in a circular 
region around each pulsar with a radius comparable to 
the LAT point-spread function (PSF) 95% containment 
as a function of energy. 

Other sources from the IFGL Catalog lfT2l either 
identified with, or associated with high-confidence to 
sources known from other wavelengths were included 
as well in the analysis model. We modeled y-ray ex- 
tended emission associated to the Cygnus Loop SNR as 
described in flTl. 

We detected significant y-ray emission associated 
with the radio shell of the y Cygni SNR and the very 
high-energy y-ray source VER J2019-H407 located on 
its rim. They were modeled using geometrical tem- 
plates, namely a 0.5° disk for the radio shell of the SNR 
lfT4l and a Gaussian for VER J2019-I-407 1 131. 

Residuals remain unaccounted for after fitting this 
model, hereafter the background model, to the LAT data 
(Fig. [TJ}). The residuals peak toward the central star- 
forming region of Cygnus X. We modeled the excess 
as a Gaussian, obtaining as best-fit parameters (/, b) - 
(79.6° + 0.3°, 1.4° ± 0.4°) for the centroid and a radius 
of 2.0° + 0.2°, with an improvement in the fit to the data 
above 1 GeV at the ~ lOcr level. The hypothesis of ex- 
tended emission provides a better fit to the data than an 
ensemble of discrete point sources. No spectral varia- 
tions across the emission region could be detected IfTOl . 



L. Tibaldo, 1. A. Grenier et al. 



3 



log(l/(W/m2/s)) -5.5 -5 -4.5 ph/bin 0.16 0.25 0.36 




Gal. longitude (deg) Gal. longitude (deg) 

Figure 1: A: 8 \im map of the Cygnus X region (MSX) showing emission from heated dust grains; massive OB stars (stars) and associations 
(circles) are outlined; the dashed circle shows the rim of the y Cygni SNR; the brightest features in the map correspond to the photon-dominated 
regions delimiting the walls of the ionized cavities carved in the surrounding dense clouds by the intense UV light. B : y-ray excess map at energies 
> 10 GeV (LAT, see text for details about the analysis). The LAT 68% point-spread function (PSF) is shown by the dashed circle in the upper 
corner. Edited from HOil . 



The residual emission (Fig. [TJ extends from the 
Cyg OB2 stellar cluster to the y Cygni SNR, peaking to- 
ward the most conspicuous clusters like Cyg OB2 itself 
and NGC 6910. However, it is more broadly distributed 
than any of these objects. It is unlikely that the broad 
y-ray excess is produced by a wind nebula powered by 
PSR J2021-H4026 (likely to be associated with the 7- 
kyr old y Cygni SNR) or by PSR J2032-H4127 (which 
is probably located behind Cygnus in the Perseus arm, 
lfT6l ). We cannot rule out the presence of a wind neb- 
ula powered by a yet undiscovered pulsar, but the mor- 
phology of the y-ray excess strongly suggests an in- 
terstellar origin. Indeed, the y-ray excess is bounded 
by the photon-dominated regions, visible through dust 
emission at 8 jam. The excess closely follows the mor- 
phology of the ionized cavities formed by the boisterous 
winds of the numerous massive star clusters in Cygnus 
X, as in a cocoon. 



3. A Cocoon of Freshly Accelerated Cosmic Rays in 
Cygnus X 

The emission from the cocoon is hard and extends 
to 100 GeV. The cocoon overlaps the TeV source 
MGRO J2031-H41 detected by Milagro. Those overlap 
as well the arcmin source TeV J2032-I-4130 detected by 
HEGRA and MAGIC. However, the flux measured by 
Milagro above 10 TeV exceeds the extrapolation of that 
measured by HEGRA at lower energies |8 1. 

Ionized gas was not part of the background model. 
It can be traced by free-free emission in the microwave 
band. Adding an ionized gas template instead of the 2° 
Gaussian provides an equally good fit, but only at the 
expense of a hard spectrum, much harder than for the 
other gas phases (§[4]). It is consistent with that obtained 
with the Gaussian template ifTOl Fig. S7]. 

In Fig. [2] we compare the y-ray spectrum of the co- 
coon with predictions based on the local CR spectrum, 
i.e. the CR spectrum in the interstellar space near the 
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Figure 2: Spectral energy distribution of the cocoon (filled circles). 
Below 1 GeV 95% confidence-level upper limits. The empty cir- 
cle gives the flux of MGRO J2031+41 subtracted by the flux of flie 
compact source TeV J2032+413() extrapolated at energies > 10 TeV 
(8]. The continuous (blue) line and upper dashed (blue) line show 
the spectrum expected from ionized gas illuminated by CRs with the 
local interstellar spectrum, assuming an effective electron density of 
f-tS = 10 cm"^ and 2 cm"^, respectively, to calculate the ionized gas 
densities from free-free emission intensities. The intermediate (red) 
dashed line shows the total IC expected if CRs had the local interstel- 
lar spectrum across the Cygnus X region, including emission from the 
tow major stellar clusters Cyg OB2 (upper -black- dotted line) and 
NGC 6910 (lower -black- dotted line), as well as from the enhanced 
optical and infrared radiation fields pervading the Cygnus X cavities 
(lower dashed line). From tlOi . 



Earth inferred from direct CR measurements, y-ray ob- 
servations of nearby clouds and radio observations of 
synchrotron emission produced by CR electrons and 
positrons. Such predictions fail to reproduce the ob- 
served y-ray spectrum. If we assume that the emission 
is purely hadronic we need to amplify the local proton 
and helium spectra by (1.6 - 1.8) X (£/10GeV)°l If 
we assume that it is purely leptonic we need an ampli- 
fication factor of ~ 60 X (£'/10GeV)°^ for the local 
CR electron and positron spectrum. The hardness of 
the y-ray spectrum points to freshly-accelerated parti- 
cles since the lifetime of TeV electrons in the enhanced 
radiation environment of the cocoon is < 20 kyr and the 
escape time for > 0.1 TeV nuclei is < 50 kyr for the 
average Galactic diffusion coefficient. 

The 7 Cygni SNR is one of the potential sites where 
the young CRs might have been accelerated. The shock- 
wave still shelters energetic particles shining in y-rays. 
CR acceleration models (taking into account CR pres- 
sure and magnetic feedback) show that, based on the 
Shockwave properties inferred from optical and X-ray 
observations, y Cygni may have released protons up to 



80-300 TeV and electrons up to 6-30 TeV at the end 
of the the free-expansion phase approximately 5 kyrs 
ago. If the dominant particle transport mechanism were 
diff'usion, and the diffusion coefficient were similar to 
that in the Galaxy at large, particles released by y Cygni 
at that time could fill the whole cocoon. Yet, there is 
no proof of a physical link between the SNR and the 
Cygnus X cavities. Moreover, the anisotropy of the par- 
ticle release with respect to the remnant (Fig.[T]i, as well 
as the comparable extension of the y-ray emission at dif- 
ferent energies, would point to an advection-dominated 
scenario instead of diff'usive transport, but there is no 
proof either that the Shockwave breaks free on its west- 
ern rim adjacent to the cocoon. In the absence of advec- 
tion the short diffusion lengths expected from the high 
magnetic turbulence level in the Cygnus X region could 
rule out the young y Cygni remnant as the only source 
of energetic particles in the cocoon. 

Superbubbles created by the collective action of SNR 
and stellar wind Shockwaves are considered as potential 
CR accelerators [e.g. J/7, J_8J. The few Myr old clus- 
ters in Cygnus X have had time to produce very few 
SNRs, if any. The total y-ray luminosity of the co- 
coon above 1 GeV, (9 ± 2) x 10^^ W at 1.4 kpc, how- 
ever, represents only ~ 0.03% and ~ 7% of the mechan- 
ical power injected by stellar winds in Cygnus OB2 and 
NGC 6910, respectively. Wind-powered turbulence in 
the high-pressure gas within the cocoon, with a typical 
energy-containing scale of 10 pc, can efficiently con- 
fine and accelerate particles. With near-saturation tur- 
bulence and a regular magnetic field of 2 nT (in pressure 
balance with the gas) the particle energy distribution 
would peak at 10-100 GeV and potentially extend up to 
150 TeV. For a typical dimension of ~ 50 pc the turbu- 
lence would scatter TeV particles for more than 100 kyr, 
in agreement with the timescale inferred from the lack 
of ^^Ni revealed by direct CR measurements IQ. 

4. The Cosmic Ray Content of the Cygnus Complex 

A y-ray excess revealed the presence of a cocoon of 
freshly-accelerated CRs over a scale of ~ 50 pc in the 
heart of Cygnus X. The average CR population of the 
whole Cygnus interstellar complex, over a scale of ~ 
400 pc, was probed by the Hi emissivity, i.e. the y-ray 
emission rate per hydrogen atom, which is derived in 
the fit of the background model to the LAT data. 

In Fig. [3] we compare the Hi emissivity spectrum of 
the Cygnus complex with the model for the local inter- 
stellar spectrum that was shown to be in good agree- 
ment with LAT measurements of diffuse y-ray emission 
from gas within 1 kpc from the Sun |fT9l . LAT data in- 
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Figure 3 : y-ray emissivity per hydrogen atom averaged over the whole 
interstellar cloud complex of Cygnus. Points include statistical uncer- 
tainties only, whereas the dashed bands show the systematic uncer- 
tainties related to the H 1 column densities and to the LAT effective 
area. The dashed (blue) curve shows the model for the local interstel- 
lar spectrum consistent with LAT observations of nearby interstellar 
clouds (I9l. From (III- 
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Figure 4: y-ray emissivity of atomic hydrogen at energies > 200 MeV 
measured by the LAT in the Cygnus region (red or light gray), the 
second Galactic quadrant 1 2 1 1 (black) and the third Galactic quadrant 
1221 (blue or gray). Results are shown for different values of the Hi 
spin temperature T$, i.e. the parameter used to determine from ra- 
dio observations the column densities of atomic hydrogen, which are 
currently the main source of uncertainty in this kind of measurement. 



dicate that the CR population averaged over the whole 
Cygnus complex is similar within ~ 20% to that in the 
much more quiet solar neighborhood (the old Galac- 
tic CR population). No counterpart to the broad ex- 
cess of y-ray emission seen at energies > 10 TeV at 
65° < / < 85° il EOl was detected by the LAT so 
far. Only the molecular clouds closer to the central star- 
forming region show a hint of a younger CR population 
from a sUght spectral hardening at energies > 10 GeV 
IfTTl Fig. 12]. This points again to an efficient confine- 
ment of the particles inside the cocoon. 

In Fig.|4]we broaden the view on the Galaxy at large 
by comparing the emissivity measurement in Cygnus 
with those in other regions of the local and outer Galaxy. 
In spite of the CR acceleration ongoing in Cygnus and 
of the steep decline in number of putative CR sources 
toward the outer Galaxy, there are no sizable variations 
of the CR densities traced by the interstellar -y-ray emis- 
sion [see also l23l . 

Unless we invoke the presence of large amounts of 
missing gas in the outer Galaxy, models assuming a uni- 
form CR diffusion throughout the Galaxy would require 
a very large propagation halo or the presence of un- 
known CR sources in the outer Galaxy to reproduce the 
LAT data Il22l . On the other hand, the strong turbulence 
associated with the stellar clusters and CR acceleration 
can increase the diffiision coefficient perpendicular to 
the Galactic disk, hence make the escape of CRs faster 
and wash out the CR density around their sources Il24l . 



5. Summary 

The Fermi LAT observations of Cygnus have added a 
few more pieces to the puzzle of Galactic CRs. The 
discovery of a cocoon of freshly-accelerated CRs in 
Cygnus X confirms the long-standing hypothesis that 
massive-star forming regions house particle accelera- 
tors. It provides an observational test case to study 
the escape of CRs from their sources and the impact 
of wind-powered turbulence on their early evolution. It 
also sheds a new light on the observation of TeV y-ray 
emission from regions harboring massive stellar clusters 
[e.g.|25]|26l|23. 

The similarity of the CR population averaged over the 
whole Cygnus complex to that in the rest of the Lo- 
cal Arm and in the outer Galaxy strengthens the dis- 
crepancy found between LAT observations and propa- 
gation models that predict a strong coupling over a scale 
of ~ 1-2 kpc between CR source densities and CR 
fluxes. Studies of Fermi y-ray observations of massive 
star-forming complexes in the inner Galaxy are required 
to further investigate this issue. 

Some of the usual assumptions concerning the injec- 
tion and transport of CRs in the interstellar space are 
challenged by the current observations. Better under- 
standing the propagation and interactions of CRs has 
important implications for the physics of the interstel- 
lar medium, hence for star formation, in the Milky Way, 
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as well as in other galaxies including starbusts. It is also 
a cracial prerequisite to look for any tracks of still un- 
known physical processes hidden in the sea of energetic 
particles filling the Universe, such as the products of de- 
cay or annihilation of dark-matter particles. All in all, 
the next century of CR Physics looks as exciting as the 
past one. 
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the Swedish Research Council and the National Space 
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